Application of density dependent parametrization models to asymmetric nuclear 

matter 
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Density dependent parametrization models of the nucleon-meson effective couplings, including 
the isovector scalar 5-field, are applied to asymmetric nuclear matter. The nuclear equation of state 
and the neutron star properties are studied in an effective Lagrangian density approach, using the 
relativistic mean field hadron theory. It is known that the introduction of a <5-meson in the constant 
coupling scheme leads to an increase of the symmetry energy at high density and so to larger neutron 
star masses, in a pure nucleon-lepton scheme. We use here a more microscopic density dependent 
model of the nucleon-meson couplings to study the properties of neutron star matter and to re- 
examine the 5-field effects in asymmetric nuclear matter. Our calculations show that, due to the 
increase of the effective S coupling at high density, with density dependent couplings the neutron 
star masses in fact can be even reduced. 
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1. Introduction 

The understanding of the properties of nuclear matter, 
at both normal and high density regions, is of crucial im- 
portance in explaining the formation and structure of 
neutron stars after the supernova explosion. The experi- 
ments with unstable nuclear beams and relativistic heavy 
ions are potential tools in determining which are the best 
equations of state (EOS) that are able to describe hot 
and dense matter. The properties of neutron stars (NS) 
are characterized by masses and radii, which are obtained 
from an appropriate EOS at high densities. The EOS 
can be derived either from relativistic or potential mod- 
els. 

The nonlinear Walecka model (NLWM) 0, Q and 
derivative scalar couplings |3j, based on the relativistic 
mean-field (RMF) approach, have been extensively used 
to study the properties of nuclear and neutron matter, 
/3-stable nuclei and then extended to the drip-line re- 
gions. In the last years some authors [1, H, H, 0, [H have 
stressed the importance of including the isovector scalar 
virtual <5(ao(980)) field in hadronic effective field theo- 
ries for asymmetric nuclear matter. The role of the 5 
meson in isospin channels appears relevant at high den- 
sities [1, 0, [H and so of great interest in nuclear 



astrophysics. 

In order to describe the medium dependence of nu- 
clear interactions, a density dependent relativistic hadron 
field (DDRH) theory has been recently suggested 0, 
11]. The density dependent meson- nucleon couplings are 
based on microscopic Dirac-Brueckner (DB) calculations 
HHEH 03 and adjusted to reproduce some nuclear mat- 
ter and finite nuclei properties [13, [LH- The main 
intention of this work is to apply different parametriza- 
tions of the density dependent meson- nucleon couplings, 
including the 5 meson, to asymmetric nuclear matter. In 
particular we will see the predictions of the density de- 
pendent coupling models when applied to the neutron 
stars (NS). In fact it is known that the introduction of 
(5-meson in the constant coupling model [8] leads to heav- 
ier neutron stars in a nucleon-lepton picture. This is not 
obvious for density dependent models. 

The paper is arranged as follows. In Sect. 2 the model 
formalism is shortly derived. The meson-nucleon cou- 
pling parametrizations are presented in Sect. 3. Results 

and discussions are given in Sects. 4, 5. 

2. The model formalism 



The Lagrangian density, with 5 mesons, 
work reads 



used in this 
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with the isoscalar (scalar, vector) er, lo p and isovector 
(scalar, vector) 5, p M , named <5, & M , effective fields. = 
dfuWv - c^a^ and G pv = - d v b p . 

The most important difference to conventional RMF 
theory is the contribution from the rearrangement self- 
energies to the DDRH baryon field equation. The 
meson- nucleon couplings g a , g^,, g p and g$ are assumed 
to be vertex functions of Lorentz-scalar bilinear forms 
of the nucleon field operators. In most applications of 
DDRH theory, these couplings are chosen as functions 
of the vector density p 2 = j p j^ with = tJjj^i/j. 

The equation of state (EOS) for nuclear matter at 
T=0 is obtained from the energy-momentum tensor. In 
a mean field approximation the energy density has the 
form [13, [n| 
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where the Fermi momentum kp i of the nucleon is related 
to its density, kp i — (3tt 2 p;) 1 / 3 . 

Since we are interested in the effects of the Nuclear 
EOS we will consider only pure nucleonic (+leptons) 
neutron star structure, i.e. without strangeness bearing 
baryons and even deconfined quarks, see [3, EH- The 
composition is determined by the requirements of charge 
neutrality and /3-equilibrium. 

The chemical equilibrium condition for a (n,p, e~) sys- 
tem can be written as 
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and the pressure is 



Ef(h) 



™2 2 

-m„o 



1 s£ 

2 m 2 



(2) 



Me = Mn - Mp = 4E sym (p)(l - 2X p ) 
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where X p is the proton fraction p p / p. The symmetry 
energy can be obtained from the energy per nucleon in 
the parabolic approximation : 



E 



sym {p) = [E/A(p, a) - E/A(p, a = 0)]/t 
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(3) 



where a asymmetry perameter a = (N — Z)/A = —p 3 /p. 
Since the electron density p e in the ultra-relativistic limit 
for non-interacting electrons can expressed as a function 
of its chemical potential, the charge neutrality condition 
is 



with E* = \jk 2 
masses are M„* 



Mr 
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i — p,n. The nucleon effective 
g<jO - gsh and M n * = M - 



g a (T + ggSs, where the scalar fields, a (isoscalar) and S3 
(isovector) are expressed in terms of the corresponding 
local scalar densities. In the pressure a rearrangement 
term appears, in the density dependent cases, as 

V R _ ( d 9a, 9o 2 , , d 9S, 96 2 
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where p 3 = p p - p n and p s3 = p sp - p sn , pi (i=p,n) and 
p s are the nucleon and the scalar densities, respectively. 

The chemical potentials for protons and neutrons can 
be written as, respectively 



p p = Jk 
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Pe = T^2 Me = Pp = X PP 



(9) 



Then, for a given p, the X p is related to the nuclear 
symmetry energy by 



3n 2 pX p 



[AE sym {p){\ - 2X p )] 3 = 0. 



(10) 



In the case of the {n,p, e~,p~) system, the con- 
stituents of neutron stars are neutrons, protons, electrons 
and muons. The threshold density for the appearance of 
muons is when the electron chemical potential is larger 
than the muon rest mass : p e > 777^=106.55 MeV. The 
chemical equilibrium for the (n,p, e~,p~) system reads 



P[i Pe — Pn Pp ■ 

The charge neutrality condition is 



(11) 



and 



Pp = Pe +P M , 



(12) 
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with the muon density p^ expressed as a function of its 
chemical potential 
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The proton fraction X p for (npe) and (npep) systems 
can be obtained by solving Eq. flOjl and Eqs. (|ll|12|13jl . 
respectively. The EOS for the /3-stable (npe) and (npefi) 
matter can be estimated by using the obtained values 
of X p . The equilibrium properties of the neutron stars 
can be finally studied by solvingTolmann-Oppenheimer- 
Volkov (TOV) equations [l|| [ljj inserting the derived 
nuclear EOS as an input. We note that the presence of 
muons slightly increases the proton fraction for a fixed 
density, making the matter softer. We will see this effect 
in the final equilibrium properties. 



Table 1. Parameters of the model 





TW [10] DDRHp 


DDRHpS 


Meson 


a 




P 


P 


S 


m, (MeV) 


550 


783 


770 


770 


980 


9i(po) 


10.73 


13.29 


3.59 


5.86 


7.59 




1.36 


1.40 


0.095 


0.095 


0.02 




0.23 


0.17 


2.17 


2.17 


3.47 


Ci 


0.41 


0.34 


0.05 


0.05 


-0.09 


di 


0.90 


0.98 


17.84 


17.84 


-9.81 



po = 0.153 fm 




FIG. 1: Density dependence of the meson- nucleon couplings. 



3. Parametrizations of the meson-nucleon couplings 

The parameters of the model include nucleon, (M = 
939MeV), and meson (m a , m u , m p , m$, see Table 1) 
masses and the density dependent meson-nucleon cou- 
plings. The density dependence parametrization used 



here, inspired by DB calculations [HI [TH, 0j| , was pro- 
posed @, 0, El as : 



9i(p) = 9i(Po)fi(x), for i = a,w,p,5, 



(14) 



with 



s l + bi(x + di) 2 . 
fi(x) = a,iexp[-bi(x - 1)] - a(x - di),i = p, (5(15) 

where x = p/ po and po is the saturation density. 
Parametrization form and parameters are taken from 
ref.[l(| for a, u mesons and from ref.@, UH for p, S 
mesons, respectively. All parameters are listed in Ta- 
ble 1. The density dependent couplings as a function of 
baryon density are displayed in Fig. [T] 
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NLp5 



FIG. 2: Density dependence of symmetry energy. Insert: pro- 
ton fraction of the corresponding /3-equilibrium (npe) system. 



For symmetric matter at saturation (pa = 0.153 /to -3 ) 
we get a binding energy E/A = e/ p — M = —16.25 MeV 
and a compressibility modulus K = 240A/eV. In order 
to remark the effects of the coupling density dependence 
we will compare the results with a non-linear (NL) rel- 
ativistic mean field model with constant couplings which 
presents very similar saturation properties (Set A of 

31.3 MeV. 
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ref . [8| ) , including a symmetry energy E, 
Both effective models, NL and DDRH are rather soft 
for symmetric matter at high density, in agreement with 
relativistic collision data, O [2(| , and Dirac-Brueckner 
expectations pll. |23[. 

As shown in refs.[J, [1] when we include the 5 coupling 
we have to increase the p coupling in order to keep the 
same symmetry term at saturation (see Table 1). Since 
at higher densities the S coupling is increasing while the 
p one is decreasing (see FigfT]), as a result in the DDRH 
choice the symmetry term will be less repulsive than in 
the NL case. This can be clearly seen in Fig. [5] at densi- 
ties above 2.5po. In the insert we present the correspond- 
ing proton fraction in a /3-equilibrated npe system. 
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FIG. 3: Equation of state for (npe) matter in different models. 
Insert: density dependence of the rearrangement terms in the 
DDRH cases 
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FIG. 4: Mass of the neutron star as a function of the radius 
of the neutron star in the two models. 



4. Neutron star results 

The /^-equilibrium nuclear matter is relevant for the 
composition of the neutron stars, as discussed in the pre- 
vious Section.. The EOS, pressure vs. density, for (npe) 
matter in the density dependent DDRH vs.NL — RMF 
models is reported in Fig. [H We see that, at variance 
with the NL results, in the DDRH cases the EOS with- 
out (5-meson is stiffer than that with the 5-meson. This 
is partially due to the softening of the symmetry term 
in the DDRHpS choice joined to a larger negative con- 
tribution to the pressure from the rearrangement term, 
see Eq.([3]), as shown inside Fig. [3l We note that both 
effects are related to the density increase of the effective 
gs coupling (see FdgfO as expected from Dirac-Brueckner 
calculations [Tl|, 0, H3 ■ 



We use the two effective nucleon-meson lagrangians, 
with and without density dependent couplings, to calcu- 
late neutron star (NS) properties, with particular atten- 
tion to the (5-field effects. The correlation between neu- 
tron star mass and radius for the /3-equilibrium (npe) and 
(npep) matter obtained by the DDRH (density depen- 
dent) and NL — i?MF(constant couplings) parametriza- 
tions are shown in Fig. 2] The obtained maximum mass, 
corresponding radius and central density for the (npe) 
and (npep) neutron star matter are reported in Table 2. 

We first note that the NLp and DDRHp results are 
rather similar, with the DDRHp interaction leading to 
a little softer matter, slightly smaller NS mass Ms and 
radius R and larger central density (see Table 2). When 
we include the 6 coupling we observe a clear effect in op- 
posite directions: the DDRH case becomes much softer 
while the NL — RMF choice shows a much stiffer be- 
havior. This can be seen from Table 2, for the variations 
in Ms/ R and central densities, but in fact it is quite im- 
pressive as it appears in Fig[4j with reference to the close 
DDRH p/NLp curves we see a clear shift to the "left" of 
the DDRHpS predictions and just the opposite to the 
"right" for the NLpS expectations. 

In general we also see, in particular from Table 2, that 
the (npe) star matter, for all models, has slightly larger 
masses and radii, and lower central densities, than the 
(npepi) star matter. This is due to the fact that the 
(npep) star matter has some larger proton fraction in the 
regions above a critical baryon density where the muon 
appears, as already noted at the end of Sect. 2. 

5. Conclusion and outlook 

All microscopic approaches of Dirac-Brueckner type to 
an effective meson-nucleon Lagrangian picture of the nu- 
clear matter are predicting a density dependence of the 
couplings. We have studied the relative effects on the nu- 
clear EOS at high baryon and isospin density, with appli- 
cation to nucleon-lepton neutron star properties. In par- 
ticular we have focussed our attention on the contribution 
of the isovector-scalar <5-meson. In fact in the "constant 
coupling" (NL — RMF) scheme the 5 leads to very re- 
pulsive symmetry energy at high density. At variance in 
the "density dependence" case (DDRH) we can have a 
"softer" dense asymmetric matter due to combined mech- 
anism of a decrease of the isovector- vector g p coupling 
and an increase of the gs (isovector scalar), which even 
leads to a larger pressure reduction from the rearrange- 
ment terms. The effect is clearly seen on equilibration 
properties of (npe) and/or (npep) neutron stars, with an 
interesting decrease of the NS mass in the DDRH case 
when the 6 contribution is included. We note that pure 
nucleon-lepton models cannot easily predict NS masses 
below two solar units. Our results seem to indicate that 
the large uncertainty of nucleon matter predictions, see 
the recent review [23[ , of relevance even for hybrid quark 
models, can be associated to the density dependence of 
the effective meson-nucleon couplings, in particular of the 
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Table 2. Maximum mass, corresponding radius and central density of the star by the different models. 





Model 


Dens. Dip 


RMF 


neutron star 


properties 


DDRHp 


DDRHpS 


NLp 


NLpS 


(npe) matter 


Ms/Mrs 


2.108 


2.01 


2.14 


2.21 




ti\km ) 


11.00 


10.29 


11.02 


11.55 




pc/po 


6.99 


7.41 


6.78 


6.44 


(npefi) matter 


Ms/Mq 


2.106 


1.98 


2.12 


2.18 




R(km) 


10.91 


10.27 


10.91 


11.30 




Pc/po 


7.14 


7.44 


6.93 


6.71 



9S- 

In conclusion we remark the interest of future work on 
two main directions: 

i) The importance of further DB confirmations of the 
high density behavior of the meson-nucleon effective cou- 
plings, in particular of some fundamental ground for the 
expected increase of the gs; 

ii) The study of dynamical effects of the isovector me- 
son fields at the high baryon and isospin densities that 
can be reached in relativistic heavy ion collisions with 
exotic beams. Differential flows and particle productions 



appear to be rather promising observables, see the recent 
refs. [13, [H, [H, H3, 111] . 
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